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Abstract 

Smoluchowski's equation is a macroscopic description of a many particle sys- 
tem with coagulation and shattering interactions. We give a microscopic model of 
the system from which we derive this equation rigorously. Provided the existence 
of a unique and sufficiently regular solution of Smoluchowski's equation, we prove 
the law of large numbers for the empirical processes. In contrast to previous deriva- 
tions we assume a moderate scaling of the particle interaction, enabling us to esti- 
mate the critical fluctuation terms by using martingale inequalities. This approach 
can be justified in the regime of high temperatures and particle densities, which is of 
special interest in astrophysical studies and where previous derivations do not apply. 
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I. Introduction 

We consider a system of dust particles of R € N different masses m\, . . . , uir, embeded in a 
(i-dimensional hot gas. Particles of size r G {1, . . . , R} are drifting according to the velocity field 
v r : R d x Rq — > R. d with a superimposed Brownian motion with diffusion constant a r E E + . Two 
particles of size r and q collide with rate a rq : R rf x Rq — > Rjj". The material coefficients e rq i G No 
determine the number of particles of size I = 1, . . . , R produced by that collision event, deciding 
for coagulation or shattering events. In this first model we take a macroscopic viewpoint, where 
two colliding particles occupy the same position in space-time (x, t), where the function a rq is 
evaluated. A complete description of the above model is given by Smoluchowski's equation [15], 
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in our case a system of reaction diffusion equations for the particle densities s r : M. d x R^ — > Rq", 
r = l,...,R (cf. [14], sect. 2) with initial conditions s° : R d -> R^: 

d t s r (x, t) = — V • (tv(x, i) s r (x, i)) + -erj! As r (x, f) 
it , i? 



-s r (x, t)^a rg (s, t) s q (x, t) + - ^ a gi (x, t) e gir s (? (x, i) i) 

9=1 ?,i=l 

s r (f,0) = s°(z) for all r = 1,... ,R (MA) 



We suppose that for any size r the particles consist of several atoms of size 1 and we set mi = 
1. Moreover, the masses of the particles are ordered as m\ < ni2 < . . . < tur. The conservation 
of the total mass M{t) = Y^r=i /r^ m ^ s r(^, i) of the system under the above dynamics, i.e. 
M(t) = M(0), t > 0, is assured by 



y~^m/ e rq i = m r + m g and a rq (x, t) = a qr (x, t) , (1) 
l=i 

for all x G M d , t G Rq", r, g G {1, . . . , R}. We also assume that e rq i is symmetric in r and q, i.e. 

The above model is commonly used to describe reaction diffusion systems and there have been 
rigorous approaches to identify equation (MA) as the limit dynamics of a suitable many particle 
system. These derivations are restricted to one-dimensional systems [1], the spatially homogeneous 
case [10] or a spatially discretized microscopic model [2, 4]. In [8] there is a derivation accounting 
for the full space dependence of the problem, using the Boltzmann-Grad limit which is applica- 
ble for very small particle densities. In [4, 10] existence and uniqueness of a solution of Smolu- 
chowski's equation are also studied. 

In this paper we give a microscopic particle model (MI) in section II, from which we rigorously de- 
rive (MA) in the spatially inhomogeneous (general) case without space discretization or restrictions 
on space dimension. Our many particle system properly describes an astrophysical system recently 
studied in [6, 14], which is explained in section V.A. It corresponds to a situation of high gas tem- 
peratures and particle densities, which is not covered by the derivation in the Boltzmann-Grad limit 
[8]. 

In this regime the dominating particle interactions are shattering collisions, so it is justified to ne- 
glect coagulation events. That means that the mass of each of the two interaction partners may 
not increase by the collision, but they are shattered into fragments of smaller or equal mass. This 
constitutes a constraint on the material coefficients e rq i given in (9), which is important to ensure 
compatibility with the microscopic particle model. Our main theorem in section III states the con- 
vergence of the empirical processes (4) to a solution of (MA) and is proved in section IV. Before 
giving a short conclusion in the last section we also discuss two apparent generalizations of the 
microscopic model (MI). 

The most important feature of our approach is the moderate scaling of the collision interaction, 
which is introduced in section II.B (M3) and discussed on a physical level in section V.A. It enables 
us to use a technique developed by K. Oelschlager [11], which was previously applied to derive the 
porous medium equation [12], or in the description of aggregation phenomena in biological popu- 
lations [9], [13]. With this technique we are able to derive Smoluchowski's equation in the spatially 
inhomogeneous form (MA), in a regime where the previous approaches cannot be applied. 



A RIGOROUS DERIVATION OF SMOLUCHOWSKI'S EQUATION 



3 



II. Microscopic particle model 

Given the macroscopic model of section I we present a corresponding microscopic many particle 
system. The most important modeling assumptions are marked by (Ml) to (M4) and are discussed 
in sections V.A and V.B. 



A. Dynamics without interaction 

Let N r (t) be the number of particles of species r G {1, . . . ,R} and N(t) = Ylr=i^r{t) the 
number of all particles at time t G M.q. The system size N is characterized by the number of atoms 
of mass mi = 1 at time t = 0: 




Let M(N, t) C N be the set of all particles and M(N, r, t) C M(N, t), r = 1, . . . , R, the subsets 
of particles of species r at time t, where each particle is identified with a unique integer number. 

(Ml) The particles are considered to be point masses with positions X^(t) G R d , k G M(N, t), 
at time t in a system of size N. Each particle of species r G {1, . . . , R} is given the rescaled 
mass mN, r = m r /N, which keeps the initial total mass Mjy(0) = J2?=i m N,rN r (0) = 1 in- 
dependent of the system size according to (2). 

(M2) Neglecting the hydrodynamic drag interaction between gas and particles, we consider the 
latter to move according to the given velocity fields v r and Brownian motion with diffusion 
constants a r , r = 1, . . . , R, introduced in the macroscopic equation (MA). 

Between two subsequent collision events the system at time t is then described by N{t) uncoupled 
stochastic differential equations: 

dX k N {t)=v r (X k N {t),t)dt + a r dB k {t) , k G M(N, r, t) , r = l,...,R (3) 

The (B k (t)) teR + , k G N, are independent Wiener processes modelling the Brownian motion of 
the particles. We always assume the existence of a filtration (!Ft)t>o, with respect to which the 
stochastic processes under consideration are adapted (cf. [11], sect. 2.B) and which fulfills the 
usual conditions [7]. 

The particle interaction is described by suitable changes of the sets M(N, r, t) and is explained in 
the next subsection. A microscopic quantity comparable to the particle density in (MA) is given by 
the measure-valued, empirical processes: 

S N , r : - M(R d ) , S N , r (t):=^ £ 6^, r = l,...,R, (4) 

keM(N,r,t) 

where A4(M. d ) denotes the space of positive, finite measures on W 1 and 5$ is the Dirac measure 
concentrated in x G M. d . Sn,t describes the time-evolution of the spatial distribution of particles 
within the subpopulation of species r. It is known by the law of large numbers that the empirical 
distribution of N independent, identically distributed random variables converges to their probabil- 
ity distribution in the limit N —> oo. In this paper we prove the convergence for stochastic processes 
which are not independent for times t > 0, due to the particle interaction. 
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B. Description of the particle interaction 

Due to (Ml) we have to specify a model for the 'collision' interaction of two point particles. 

(M3) We take a stochastic model determined by a rate depending on the distance of the interaction 
partners k and I. The scaling of this rate is given by 

^W N (X k N (t) - X l N {t)) , whereW N (x) = a d N W 1 (a N x) , (5) 

with aw = N@/ d and a moderate scaling parameter < (3 < 1. We assume that W\ is sym- 
metric and positive with ||Wi||i = 1. It follows that ||W}v||i = 1 for all N G N and 
limjv->ooWjv(- — x) = 5% for all x G R d in the sense of distributions. 

In contrast to the usual hydrodynamic scaling with (3=1 this leads to a microscopically large 
interaction volume. This assumption is motivated and justified in a physical context in section V.A. 

(M4) Instead of considering pair interactions (see sect. V.B) we assume that every particle k G 
M(N, r, t) interacts with an effective field of all other particles of species q with rate 

a N>rq {X k N (t ) , t ) : = min { C a , a rq {X k N (t ) , t ) 

((SW,,(t) * W N )(X k N (t))-5 r , q W N (d)/N^}, (6) 
where the macroscopic collision rate given in section I. 
In (6) we used the generalized convolution product 

S N , q (t) *W N := [ W N (. - x) S N , q (t)(d d x) = 1 V W N (. - X l N (t)). (7) 

jRd ' N leM(N, q ,t) 

By substraction of the term including Kronecker's delta in (6) self-interaction is excluded. The 
rate is bounded uniformly in N by a suitable constant C a , which is specified in condition (C5) in 
section III.B. This cut-off prevents diverging interaction rates due to high particle concentrations in 
the limit N —>■ oo. Each possible interaction event is described by a jump process 

: = PN,rq XM{N,r,s){k) CL N , rq (X%( S ), s) ds^j G {0, 1}, (8) 

where 0^ rq : M^" — > No are independent standard Poisson processes with a transformed time argu- 
ment in the brackets (cf. [11]) and \A G {0, 1} is the indicator function of the set A. 
The process a*^ k rq jumps from to 1 at some time t > if particle k exists in t— , belongs to species 
r and interacts with a particle of species q at time t. After the interaction the number k is removed 
from the sets M(N, r, t) and M(N, t). The mass of particle k is distributed on the interaction prod- 
ucts according to the microscopic material coefficient e rq i G No- The latter fulfills conservation of 
mass and is related to its macroscopic counterpart e rq i in the following way: 

R 

^m/e r(? / = m r , e rq i = e qr i = e rq i + e qr i for all r, q, I = 1, . . . , R. (9) 
l=i 

We note that this also constitutes a condition on e rq i, corresponding to the absence of coagulation 
mentioned in section I. The particles resulting from the interaction are located at X^(t) and obtain 
new numbers starting with maxjj? G M(N, s) : s < t} + 1, which were previously not assigned 
to any particle. These numbers are added to M(N, t) and the subsets corresponding to the various 
species. We note that any process a*^ h rq only jumps once, since after that jump the respective particle 
k disappears, i.e. XM(N,r,t+.)(k) = 0. 
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C. Complete description of the model 

Using a generalized L 2 -scalar product we can formulate the time evolution of the empirical 
processes in a weak sense. For all / € C 2 (K. d , R) and r = 1, . . . , R we have 



(S N , r (t)J):= [ f(x)S N At)(d d x) = ± Yl /(*&(*)) = 

jRd keM(N,r,t) 

= jr E f(x k N m+ fdl^ Yl /(*&(*))]■ do) 

keM(N,r,0) \ keM(N,r,s) J 



Inserting the expression for dX^(s) from equation (3) and using Ito's formula [7] we get: 

(%r(t),/) = ^ £ f(X k N (0)) + ^ f Y Vf(X k N (s))-dB k (s) 

keM(N,r,0) keM(N,r,s) 

+^ f E • -) + f A/(^( S ))) 

- 70 keM(N,r,s) V ' 

-^eT e 

9=1 • /0 keM(N,r,s) 

+^E / E w^(*o ai) 

g ,/=l- /0 keM(N,q,s) 

The first integral term describes the stochastic fluctuations of the particle positions and the second 
one particle transport and diffusion, resulting from the interaction free description (3). The next 
two terms consider the change of the sets M(N, r, s) in (10) due to the loss of particles of species r 
after interactions with others, and the gain of such particles from products of other interactions. We 
separate the fluctuation terms due to stochasticity in the free particle dynamics and the interaction 
in stochastic integrals. So we get for all / € C 2 (R d , R) and r = 1, . . . , R the complete description 
of our microscopic model: 

(S N , r (t), f) = (S N , r (0), f) + jf (s Ntr (8), V/ ■ V r (; *) + yA/) ds 

R i-t R pt 

~E / ( S NA S )'f a N,rq(;S))ds+y2 / , / e qlr a N ,ql(- , s)) ds 

q=l J ° q,l=l J ° 

+M 1 N , r (f, t) + M 2 N * r (f, t) + M 2 N b r (f, t) , With 
MkAM = Jjf E Vf(X k N (s)).dB k (s) 

J ° keM(N,r,s) 

M 2 N a Af,t) = ~Yl I E /(^(s))(«^(^)-a^(^( S ), S )^) 

q=l J ° keM(N,r,s) 

(W = Jf E / E f(X k As))e qlr (a% k gl (ds) -a N , ql (X k N (s),s)ds) 

q,l=l J ° keM(N,q,s) 



M 2b 
1V1 N 
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and initial conditions <Sjv,r(0) = N 1 J2keM(N,r,o) S x^(o)- ( Mr> 

This set of equations combines all features mentioned in the preceding two subsections and is 
used to derive the macroscopic model (MA), shown in the next section. 



III. Derivation of Smoluchowski's equation 

We show how to obtain (MA) heuristically from our microscopic particle model (MI), leading 
us to a proper formulation of the main theorem. 



A. Heuristic derivation of the macroscopic equation 

The empirical processes Sn jT are defined as solutions of (MI). For this subsection we assume 
that for every r = 1, . . . , R they converge to limit processes S r : [0, T] — > A4(R rf ) on a compact 
time interval [0, T] in a yet unspecified sense. The limit processes are assumed to be absolutely 
continuous with respect to Lebesgue measure on R d and therefore have densities s r , which should 
be in C 2 (R rf x Rq",Rq~). With the generalized scalar product defined in (10) we therefore have 
(S r (t),g(.,t)) = (s r (.,t),g(.,t)) for all g G C b (R d x [0,T],R). We also assume the validity of 
conditions (CI) to (C7) given in the next subsection. 

In section IV.D we get the following for the stochastic integrals in (MI) for any T > 0: 



lim E 

iV->oo 



sup \M^ r (f,t)\ 

te[0,T] 



= fori = 1,2a, 26 , (12) 



so the fluctuation terms asymptotically vanish in any compact time interval and the limit equation is 
supposed to be deterministic (see (14)). The convergence of the Sjv, r should be sufficiently strong 
to assure the following: 

lim E[(S N , r (t),g(.,t))] = (s r (.,t),g{.,t)) 

lim E [(S N , r (t), {S N , q (t) * W N ) g{.,t)}} = (s r (.,t), s q (., t) g(.,t)}, (13) 

for all r, q = 1, . . . , R, t G [0, T] and g G C b (R d x [0, T], R d ). The first condition assures the con- 
vergence of the drift and diffusion term in (MI) and the second one is needed for the interaction 
terms. We formally substitute the above limits into (MI) and notice that the self interaction term in 
(6) vanishes for — > oo. Therefore we get the following deterministic integral equation for all test 
functions / G C 5 2 (R d , R), t G [0, T] and r = 1, . . . , R: 



(s r (.,t), f) = (s r (.,0),f) + f^ds (s r {., s), V/ • v r (., s) + ^ 2 a/ 

R rt 

-V] / ds(s r (.,s),fa rq (.,s)s q (.,s)) 

q =l J ° 
R rt 

+ V] / ds{s q (.,s),fe q i r a q i(.,s)si(.,s)) (14) 

q,l=l 

After partial integration in the transport and diffusion terms one immediately recognizes this as a 
weak version of Smoluchowski's equation. Using (1) and (9) it is easy to get the last line in the 
form (MA). 
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Therefore we showed that, assuming the empirical processes converge, their limit densities fulfill 
a weak form of Smoluchowski's equation. In the next subsection we explain how to prove this 
convergence in an appropriate rigorous limit sense, which can be seen from (13) to be of L 2 -type. 

B. Convergence theorem 

To formulate the convergence theorem we use the following distance function between the em- 
pirical processes (MI) and the solution of Smoluchowski's equation (MA) specified in (C3) below: 

dN,r(-, t) := h N ,r(-,t) - s r (., t) , where h N , r (.,t) := S NjT {t) * W N , (15) 

for all r = 1, . . . , R and t G Rq". The convolution kernel Wn smooths out the empirical processes 
and obeys the following regularity conditions: 

(CI) Wn is a different scaling of the interaction function Wn and both have to fulfill: 

Wn(x) = a N Wi(aNx) and Wn(x) = a N W\{a.Nx) , where 

&N = N m^ < p < d and aN = N p /d ^ <(3< P_ 

d+2 d+l 

The scaling parameter [3 plays no role in the dynamics of the many-particle system. However, by 
the above assumptions some restrictions on the parameter j3 determining the moderate interaction 
are introduced. 

(C2) The unsealed function W\ € L l n C 2 (IR d ,M) is symmetric, positive and standardized, i.e. 
|| Wi||i = 1- We also need J Rd \x\ W\(x)d d x < oo and the Fourier transform W\ has to fulfill: 

a) W x G C b 2 (R d ) 

b) |Wi(r)| < Cexp(-C"|f|) 

c) |AWi(t)| < C(l + |r| 2 )|iyi(f)| 

d) v |Wi(t;r)|, v > 0, monotonicly decreasing for all fixed f € M. d 

A Gaussian probability density is an example for W\ which obeys these conditions. To the knowl- 
edge of the authors there is no proof of the existence of a sufficiently smooth solution of the macro- 
scopic equations, therefore we have to assume the following: 

(C3) There exists a positive, unique C 2 -solution (s\, . . . , sr) of Smoluchowski's equation (MA) 
in the time interval [0, T*] for some positive T*. The functions s r (.,t) and their partial deriva- 
tives are L 2 (R d , M+)-bounded uniformly in t G [0, T*]. 

(C4) The macroscopic collision rate d rq {.,t) given in section I should be Lipschitz continuous, 
bounded and fulfill the conditions (1) for all t G [0, T*]. The macroscopic material coefficient 
e rq i should obey condition (1) and together with its microscopic counterpart e rq \ given in 
section II.B, it should fulfill (9) and be symmetric in r and q. 

(C5) The upper bound C a for the microscopic interaction rates (6) is given so that the limit equa- 
tion is not affected, C a > max r ?e {i ) sup te [ 0jT .] ||a r . g (.,t)s g (.,t)|| 00 . 

(C6) The velocity fields of the different particle species have to fulfill 

v r G Cl(R d x [0, T*],R d ) for all r = 1, . . . , R. 

(C7) The diffusion constants of all particle species have to be positive, i.e. a r > for all r = 
1,...,R. 
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We note that our proof only applies if all particles are Brownian. Now we are ready to formulate 
our main convergence result. 

Ef=ill^(.,o)||i 



Theorem. With conditions (CI) to ( C7) and Xim^^^E 
R 



it is 



lim E 



V SUp \\d Nyr (.,t)\\l 
T =l *G[0,T*] 



0. 



(16) 



Convergence at time t = is given if the initial conditions of (MA) and (MI) are compatible. One 
possibility is to take the particle positions X^(0), k <G M(N, r, 0) as independent, identically dis- 
tributed random variables with suitably normalized densities Sr/X^-Li ( s r> 1) f° r all r = 1, . . . , i?. 
For discussion of this point see [12] (sect. 4B). 

To formulate the result without the smoothing convolution kernel Wn we introduce a metric on 

M(M. d ) by 

D(p,u) := sup{\{p -v,f)\:f€ tt D } for all p, v e M(M d ) and 

n D := {/ G C l b n L 2 (M. d ,R) : ||/||oo + ||V/||oo + ||/|| 2 < 1}. (17) 

This quantifies a distance between the empirical processes defined in (MI) and the processes S r (t) := 
J s r (x, t)d d x given by the solution (C3) of the macroscopic equation. As the theorem states con- 
vergence in an L 2 -sense the convergence in the weak sense (17) is easy to conclude. 



Corollary. With the conditions of the theorem we have 



lim E 




sup D(S N>r (t),S r (.,t)) 

t€[Q,T*] 



0. 



(18) 



IV. Proof of the convergence result 

A. Preliminaries 

The following lemma is useful in central estimates of section IV. 
Lemma. With f e L 2 n C£(R d , M), V/ G L 2 (R d , R d ) we have 

||/-/*^||2<C« W 2 ||V/||^ and \\f -f*W N \\ 00 <Ca N 1 \\Vf\\ 00 . (19) 

An analogous estimate is true, if Wn and a m are replaced by Wn and cxn- 
For any finite, positive measure fi on M. d and with Un(x) := \x\ Wn(x) there is 

\\fi*U N \\l <Ca 2 ^- 2 \\n*W N \\ 2 2 + (ijl, l) 2 exp(-C'a e N ) foralle>0. (20) 
For any finite, signed measure p, on W 1 it is 

\\n*W N \\l<\\n*W N \\l (21) 
Proof, see [11] sect. 4A,B and [12] sect. 5B, or [3], sect. 4.3 

In the proof of the lemma there is essentially made use of the conditions (CI) and (C2) on the 
interaction function Wn and the kernel Wjy. Due to the conservation of mass in the microscopic 
system (9) and with (2) we get the following bound on the empirical processes, 

R 1 N 

r=l keM(N,t) 
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We also use the following property without explicitly noting it for all suitable / and g, such that the 
expressions are well defined: 

(f,g*W N ) = (f*W N ,g) and (/, g * W N ) = (/ * W N , g), (23) 

because W\ is symmetric according to (C2). Throughout this chapter C, C etc. denote suitably 
chosen constants, whose value can vary from line to line. 



B. Proof of the theorem 

To prove statement (16) we first look at the time evolution of the quantity 

\\dN,r(-Ml = ||^JV,r(-,*)||| - 2</l JV)? .(.,t), J 9 T .(.,t)> + ||s r (-, *) |||- (24) 

The dynamics of the first two terms is obtained analogous to (11) using (3), (8), (23) and Ito's 
formula: 

\\hN,r(;t)\\ 2 2 = J^ E * W N )(X k N (t) - X l N (t)) 

k,leM(N,r,t) 

(h N , r (.,t),S r (.,t)) = 1 Yl M'Q * W N )(X%(t)) 

keM(N,r,t) 

We just have to replace the test function / in (11) by £jv,r(*) * (Wn * Wjv) res P- s r (., t) * Wn- 
The expansion of the third term in (24) follows from the macroscopic equation (MA): 

\\s r (.,t)\\ 2 2 = K(.,0)||! + J ds(s r (.,s),-2V ■ (v r {.,s)s r {.,s))+a?As r {.,s)} 

+ V / ds ( s r (., s), -2s r (., s) a rq (., s) s q (., s) + s q (., s)Ve, lr a ql (., s) s t (., s) ) 
q =i Jo \ l=i I 

Combining the parts suitably by using (9) to express e q i r in terms of e q i r we get: 

H<*iV,r(.,*)ll2 = l|rfjV,r(-, 0)||1 

+al \ ds (d Nyr (.,s),Ad Nyr (.,s)) 
Jo 

^(S N , r (s),V(d N , r (.,s) * W N ) ■ v r {.,s)) - (s r (. , s) ,V d N , r (., s) ■ v r (.,s))^j 



+2 [ r/.s 
Jo 

R rt 



-2V / ds ((SW,r(s), (d N>r (.,s) * W N )a N}rq (.,s)) 
-{s r (., s), djv,r(-, s) a rq (., s) s q (., s)}) 

R ft ( 

+2 V ds e q i r [{S Ntq (s), (d N ,r{-,s) * W N ) a N>q i(., s)) 

-{s q (., s), d N>r (., s) a q i(., s) si(., s)}) 
+2 [^ E V(d N A,s)*W N ){X%(s)).dB k {s) 



keM(N,r,s) 
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R rt 



" 2 E / E (dAr,r(.,a)*^)(^( S ))(a^,(cfa)- 0^,(^(8), s) da) 

9=1 ^° keM(N,r,s) 

+2 E / ^ E (^, r (.,s)*^iv)(^( S )) e ^(a^(d S )-a^(^(5)^)^) 

oJ=r° keM(N,a,s) 



q,l=l u k&M{N,q,s) 

■ / (is (%r(s), 1) 
JO 



+- — 


iV 


! (W N 


* Wjv)(0) 




iV 




*WV)(0) 




N 2 




* w N )(6) 



N 

/ J I ds (S Ntr (s),a Ntrq (.,s)) 

a=l J ° 



R r t 



\~] / ds e q i r {S Nyq (s),a N ^ q i(.,s)) 
q ,l=i Jo 

R rt 

E/ E \ a *N k r q ( ds ) -a N , rq (X N (s),s)ds) 

(7=1 keM(N.r.s) 



= \\d N ,r(-> °)Hi + J (T^ r (s) + T^ r (s) + T^ r (s)+T^ r (s)) ds 
+2 (i\4 >r (t)+M^ r (t) + M» r (t)) 

+T° ir (t) + 7$ r (t) + T$ r (t) + A4$ r (t) + M°» (25) 

The terms in the above sum are labeled line by line. Tjy derives from the diffusion due to Brownian 
motion, T N from the particle transport, T N from the loss and T N from the gain of particles due 
to interacions. The stochastic integrals M^ r , M N % r and Mff r represent the fluctuations due to 
stochasticity in the free particle dynamics and the interaction. The remaining terms are corrections 
resulting from the expansion of ||/ijv, r (., With Ito's formula and (3) we get for T Nr : 

^ Y, d(W N *W N )(X N (t)-X l N (t)) = 

k,leM(N,r,t) 



N 2 

k,l<=M(N,r,t) 

+ 



^A(W N * W N )(X k N (t)-X l N {t)) (dt + dtfj 



= 2{S N , r (t), V(h N , r (.,t) * W N ) ■ V r {.,t)) dt + ^ ^ V(h N , r (.,t) * Wjv) • d& (t) 

keM(N,r,t) 

+a 2 r (h N:r (.,t),Ah Ntr (.,t)}dt-^ Yl A(W N *W N )(6)dt 

keM(N,r,t) 

T^ a r and T N b r can be derived analogously by considering the change of the sets M(N, q, t), q = 
1,...,R. The fluctuations of these corrections are separated in stochastic integrals M^^, and M^ r - 
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Estimate of the correction terms: 
With (CI), (C2) and (22) we get: 

K, r (t)\ < ^ ( [ \W N {x) AW N (0 - x)\d d x) sup(S N , r (s), 1) < 

iV \JR d J s<t 

<° t [ N^ 1+2 ^\\AW 1 \\ oc \W N (x)\d d x < CtNfo+W- 1 

This term vanishes in the limit N — > oo, because with (CI) it is $ < dj (d + 2). Using also (6) and 
(9) we get completely analogous: 

\T° N - r (t)\ + \T^ r (t)\<CtN^ 
The stochastic integrals are estimated in section IV. D. 

Estimate of T^ r (s) = a 2 {d N , r (., s), Ad N>r (., s)) 

After partial integration we get with condition (C7): T^ r (s) = — a 2 \\VdN,r(-, < 0- This 
term is negative and can be used to cancel positive contributions of the same kind arising in the 
estimates of r and r . 

Estimate of 

T ^r( S )= 2 ((%r(s),V(d7V,r(., s)*W N )-V r (., s)) ~ (s r (., S ),Vd N , r (., S )-V r (., s)>) 

To contract the two brackets we make the following replacement: 

{S N , r (s), (VeZjv, r (-, s) * W N ) ■ v r {., s)) = (h N , r (., s), Vd N . r (., s) ■ v r (., s)} + R 2 Ntr (s), 
where the correction term is estimated using (C6): 



\R%Js) 



(^S N ,r(s), J ^ d d uW N (u)Vd Njr (. -U,s) ■ (V r (-,S) - V r (. ~ U, s))^ 



< 



loo < 



< /5Ar ) r(s),j^ d ^«Wj V (^|u||Vdjv,r(.-«,s)|^ ||W r (.,, 

< C(% r ( S ) * U N , \Vd N , r (.,s)\) < C (C\\S N , T {S) * U N \\ 2 2 + i\\Vd N , r {.,s)\\i\ , 

with Un(x) := |x|Vl^v(x). This is true for all C > using Cauchy's inequality. With the second 
statement of the lemma (20) and (22) we have for all e > 

\R 2 N , r ( s )\ < cc-^^ii^c-.s)!!! + <5 JV>r ( a ),i> 2 e- C7 ' a ^) + ^Hvdjv.rC-,*)!!! < 

< CC(c%- 2 {\\d N , r {.,s)f 2 + C) + e- c "^) + -^\\Vd N:r (.,s)\\l 
using (C3), because with the triangle inequality it is 

< (\\d N , r (;S)h + \\Sr(.,s)hf < 211^0,5)111 +2||s r (.,s)||i. 

C is chosen after the estimate of Mh , so that the term arising there and Or- 1| Vdjv> (•, s) \\\ cancels 
with the negative contribution from the estimate of T^ r . Choosing e = \ the constant terms in the 
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above estimate of \R% r {s)\ and the prefactor of ||djv,r(-> s)||| vanish in the limit N — > oo. Now we 
can write 

T^ r (s) = 2(d N , r (; S), (Vd N , r {., S)) ■ V r (., s)) + 2i& >r (s). 

With (C6) and the estimate 

(d N)r (.,s), (Vd Njr (.,s)) ■ V r (.,s)) =- (W^ r (., s), V r (., s)) 



; (di r (.,s),V-tv(.,s)) 



< TTll V • *V(-,s)||oo||diV,r-(-,s) 



we get after a suitable arrangement of constants and terms: 

|T£»| < C(l + C)\\d N , r (., s)\\t + Sr\\Vd N , r (., s)f 2 + C 0(N-^ d ), 

because the iV-dependent prefactors of ||<iiv,r (•? s) ||| vanish monotonically with N —> oo and N~@/ d 
is the leading order in iV of all constant terms. 

Estimate of T% r {s) = -2J2f =1 ({S N;r (s), (d N>r (., s) * W N ) a N>rq {., s)) 

-{s r (., s),d N;r (., s) a rq (., s) s q (., s))) 

First we make the same substitution as before: 

(S N , r (s), (d N!r (., s) * W N ) a N;rq (., s)) = (/ijV,r-(-, s),d N;r (., s) CL N , rq {., s)) + R 3 ^ rq (s), 

and we can estimate the correction term analogously to R% r using (6), (CI), (C2), (C4), (20) and 
(22): 



S N , r (s), / d d uW N {u) d N , r {. - u,s)(a N: r q (.,s) - a N , rq (. - u,s)) 

<(S N ,r(s), / d d uW N (u)\u\\d N:r {. -U,S)\ 
\ JR d 

(||a r? (,s)|| 00 ||V(5 A r, g (s)*^)|| 00 + ||Va r? (.,s)|| 00 ||5 A r, g (s)*^|| 00 ) < 



< 



< (SnAs) * U N , |djv,r(.,*)|> C (af 1 + a^XU^IU + < 



< Ca 2d+2 (a 2 ^\\\d N A.,s)\\l + C) + e- c "^) + \\d N , r (.,. 

For the first term in the above estimate to vanish in the limit N — ► oo, we choose e > so that 7 := 
0(2d + 2) + (3(2e - 2))/d < 0. That means < e < 10 - (3(d + 1)), which is possible due t 

condition (CI) on the scaling parameters (3 and (3. Now we look at the remaining term in T^ r : 

R 

-2^Z (i h N,r(-, s),d N , r (; S) a N , rq (., S)) - {s r (., s), d N , r (., s) d rq (., s) S q (.,s)fj 



q=l 



To contract the two brackets we have to compare the microscopic and macroscopic interaction rates: 

(hN,r{-, s),dN,r(-, s) CLN,rq(-, s)} ~ (s r (., s) , djV,r(-, s) & rq (., s) S q (., s)) = 



&N,rq 

(., S ))+R% rq {s), 



with the correction term 

3 3f> 



R N,rq( S ) = ( S A; s),d N ,r{-, s)(a N>rq (., s) - d rq (., s) S q {., s))) . 
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With the definition of ajv jrg (., s) in (6) and (C5) we get the following estimate: 

\R 3 N,rq(s)\ < ^S r .(.,s),|d A r, r (., S )||a rg (.,s)| ^^(.,5) - (s) * T^V | + ^ T^V (0) 

■ S N>q (s) * W N \\ 2 



< 



s 2 * 



< \\s r (-, s) a rq (-, s)\\ O0 \\d N) r(., 

+5^-^(0) {\\a rq (., s)\\ 2 2 + \\d N;r (., s)\\l) 

It is \\s q (.,s) - S N) g(s) * W N \\ 2 < \\(Sq(.,s) - S Nj q(s)) * W N \\ 2 + \\s q (.,s) - S q (.,s) * 

and with the third statement (21) of the lemma 

s) - S N , q (s)) * W N \\ 2 < \\(s q (., s) - S N , q (s)) * W N \\ 2 < 

< \\dN,q(;S)\\ 2 + ||Sq(.,s) - S q (.,s) * W N \\ 2 . 

Therefore we get with the first statement (19) and (C3) 

\R%rq{s)\ < C(\\d N , r {;S)\\l + H«W->*)ll2 + C (aj* + &~ N 2 )\\V S q (., S )\\ 2 ) 

+5 r>q C"NP- 1 (l + \\d N A,s)\\ 2 2 ). 
After treating the correction terms we get for the main contribution using (6): 

\(d N!r (.,s),d Njr (.,s)a N ,rq(;s)}\ < C a \\d N , r {., s) ||| 
Arranging all terms analogously to r we finally have the estimate 

R 

\T%,r( 8 )\ <Cj2W d ^ S )\\2+°( Nl ) . With 7 <0. 



W N \ 



Estimate of T^ r (s) = 2 £ ? ,;=i e gir f (S Njq (s), (d N , r (., s) *W N )a Njql (.,s)) 



■(s q (., s),d Ntr (., s) a q i(., s) si(., s))j 



Obviously, the estimate of this term is completely analogous to the one of with the same 
result except for different constants, 

R 

\T^ r (s)\<Cj2\\dN, q (^)\\ 2 2 + 0(N^ , with 7 <0. 

9=1 

Putting all the estimates together, taking the absolute value, the supremum over all t € [0, T] for 
some T < T* (see (C3)) and the sum Ylr=i on b otn s ^ es °f equation (25) we arrive at: 

R / 



sup \\d N , r (.,t)\\ 2 2+(4-^] I ds\\Vd N , r (., 
1 te[0,T] V CJ Jo 



r=l 



R 



<V [\\d N A;0)\\ 2 2 + TC SU P \\d Ntr (.,t)\\l] + (1 + C + T)0(W) 

' te[o,T] J 

R 

+^ sup (2|ikfAr,.(*)l + 217^^(01 + 2|ikT^(t)| + + I^VC^l) 

^e[o,T] v > 



(26) 



The choice of C below ensures the positivity of all occuring terms. The leading order of all constant 
terms, including the estimated correction terms, is characterized by 7 < defined in the estimate 
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of T]y r . Taking the expectation on both sides of (26) we can use the estimates (35) and (36) of the 
stochastic integrals in section IV. D and get: 



E 



R 



SUp ^(.Ml+ivl-CN^ 2 -Z\ f ' ds\\Vd N , r (.,8)\\\ 

1 te[0,T] V C J Jo 



r=l 



< 



< E 



R 



J2(\\d N A-M\ 2 2+TC" sup {\\d N , r (.,t)\\t) 
~~ 1 te[o,T] 



+ (1 + C + T) O (N^) , (27) 



where the leading order of constant terms remains unchanged with (CI). Now we can choose C, oc- 
curing in the estimate of T^ r , and Nq G N large enough, so that the pref actor of 

Jq ds ||Vdjv,r(-) s )||| on th e lefthand side of (27) is positive for all N > No (see (C7)). Conse- 
quently this term can be neglected and after a short rearrangement we have for all N > Nq and 
< T < 1/C", where C" is taken from (27): 



E 



R 



r=l t6[0,n 



< 



l-TC" 



E 



+ (l + T) 0(J\P 



(28) 



Taking the limit N — > oo on both sides, the constant terms vanish as 7 < and convergence at t = 
is given in the theorem. This finally proves statement (16) for all t € [0, T], but the above constraint 
on the time interval is not essential. At t = T all the conditions for the theorem are fulfilled as long 
as T < T*, which enables us to apply the proof again with the same constants arising. The length 
of the next time interval is subject to the same constraint and after a finite number of repetitions this 
proves the theorem. 



C. Proof of the corollary 



We have the following estimate: 

\(S N , r (t)-S r (.,t),f)\ < \(S N , r (t)-h N , r (.,t),f)\ + \(d N , r (.,t),f)\ < 

< \(S N , r (t), f-f*W N )\ + \\f\\ 2 \\d N , r (.,t)\\2 < 

< (S N , r (t), 1)11/ - / * ^vlloo + |m| 2 ||aM.,i)|| 2 

Therefore we get for all / G Op (see 17), using statement (19) of the lemma and (22): |(<SV,r(i) — 
s r (.,t),f)\ < CdlV/Hooci^ 1 + ||eZjv,r-(-,i)ll2)- Hence we have for the metric D 



\D(S N , r (t),S r (t))\ < Ci&rf + WdNri.rfh) for all r = 1,...,R. 
It is straightforward to see that the statement (16) of the theorem implies 



(29) 



lim E 



V SUp \\d Ntr (., 
r=l *e[0,T*] 



0. 



Using this, equation (29) immediately gives 

R / \ 
E sup D(S N , r (t),S r (t)) 



lim E 



< C lim E 



R 



V SUp ||djv,r(-,*)l|2 

r=1 te[o,T*] 



0, 



which proves the statement (18) of the corollary. 
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D. Estimate of the stochastic integrals 



In the following we estimate the stochastic integrals occuring in (MI) and (25). First we show 
for every fixed system size N € N that the integrands and integrators fulfill the necessary regularity 
conditions so that standard techniques of stochastic integration taken from [5, 7] can be applied. 
This part is kept short and can be read in more detail in [3] (sect. 3.1.1, appendices B and C), 
following the work in [11]. According to those results we find an estimate uniformly in N using 
Doob's inequality [7]. First we consider the terms 



N 



oo d 



EE / 'Wm^w 
k=i i=i j ° 

where I l N ^ kii {f\t) := XM(N,r,t)( k ) 

R oo „t 

M 2 N a , r (f,t) = "^EE J T N,r,kU,s) dP^ rq (8) 



q=l k=l ' 

where lft, >k (f,t) := X M(N,r,t)(k) f(X k N (t)) 

and P^ rq (t) := af (t) - [ X M(N,r,s)( k ) a N ,r q (X^(s), s) ds (30) 

Jo 

occuring in (MI). Mffl r can be handled analogously to Mff r . The integrators B\ and P% rq in (30) 
are Brownian motions and time-inhomogeneous compensated Poisson processes (8) with bounded 
rates (6). Therefore they are square integrable martingales with respect to the filtration (Ft)t>o> 
mentioned after equation (3). The integrands l\ rki and ljf r k can be replaced by their left contin- 
uous versions as we have 



E 



/V^(/^)-4,^(/^-)i 2 ^)(«) 

Jo 

l\ll\ k {f,s)-ll\ k {f,s-)\ 2 d{P^ rq ){s) 
Jo 



0, 



because the quadratic variational processes (Bf)(s) = s and (-PjvrgH 5 ) = Jo XM(N,r,u)(k) 
aN,rq{X%{u),u) du are continuous in s. On Poisson processes with time dependent rates see e.g. 
[5], theorem II.3.1 on page 60. Therefore we can assume the integrands to be predictable processes 
with respect to (Ft)t>o- For every fixed system size N G N it is 



E 



l4r A i(/.')W)W 



E 



l N,r,k 



(f,s)\ 2 d{P*)(s 



< OO, 



(31) 



using / £ Cl 



and estimates analogous to (32) and (33) below. Hence M^ r (f,t) and 



M^rif, t) are martingales with respect to (Tt)t>o- On the other hand we have for all r 
t > and N € N: 

E[\M l N , r {f,t)\ 2 ] = 



, R, 



E 



W E 



k€M(N,r,s) 
rt d 



leM(N,r,s) 



£ £ d x J(X%{s)) d Xj f(X l N (s)) d(BtB]){s) 



k,l€M(N,r,s) i,j=l 
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W E 



\Vf(X k N ( S ))\ 2 ds 

keM(N,r,s) 
t 



(J 



2 ft 



^ / E[{S N , r (s),\Vf\ 2 )}ds< 



<^I|V/||^ / E[{S N As),l)]ds<C^\\Vf\\l t<™ 



(32) 



E[\M 2 N * r (f,t)\ 2 ] = 



1 R 



N 2 

p,q=l 
1 R 



f(X k N (s))f(X l N (s)) d(P^ rp , P l N , rq ) 



=i 



<- 



k,leM(N,r,s) 

£ f 2 (X k N (s))a N , rq (X k N (s),s)ds 

keM(N,r,s) 

SUp y^||aAT,rq(-,s)||oo E [{S N>r (s) , 1)] ds 
se[0,t] j / Jo 



< 



< 



N 



J < oo 



(33) 



This is true with (22), / G C 2 (R d , R), (C2) and (C4). After an estimate of M 2 ^ analogous to M 2 ^ 
we can apply Doob's inequality and get for i = 1, 2a, 26, r = 1, . . . , R and all T > 0: 







2 






sup |M^ r (/,t)| 


< # 


sup |M^ r (/,t)| 2 




te[o,T] 




te[0,T] 



<4E[\Mh tr (f,T)f 



(34) 



so the result (12) follows from (32) to (34). 

For the terms occuring in (25) we have M l N r (t) = M l N r (dN, r (-, t) * Wat, t) for i = 1, 2a, 26 and 
therefore a similar reasoning applies in this case. To get an estimate in the limit N — > oo we note 
that with (CI) and (C2) Wn(x — .) is a probability density on M. d for all x G R d . Hence we have 
for all / : R d — > R, where E 1 ^ denotes the corresponding expectation value: 



S \|2 



\(f*W N )(x)\ 



[ W N {x-y)f{y)d d y 



E WN [f] 2 <E w Jf 2 ] = (\f\ 2 *W N )(x) 



J w N 



We also have with (CI), (C2) and (22) for all t G [0, T*] and r = 1, . . . , R: 

||^,r(.,t)||oo < E 11^1 («JV(- -^(*)))IU < CiY^% r (t),l) < CN? 



N 

keM(N,r,t) 

With the last two considerations we get for all T G [0, T*], starting with the third line of (32) and 
second line of (33) and using (23), 



E 



<^E 
~ N 



(S N , r {s),\V(d N , r (., S )*W N )\ 2 )ds 



< 



< ^E j\h N , r (., S ),\Vd N , r (.,s)\ 2 )ds 



< CN^olE 



\\Vd N , r {.,s)\\ 2 2 ds 



uo 



E 



R 



\M% a AT)\ 2 <^ E 



{Sn,t{s), \dN,r{-, s) * Wn\ 2 a,N,rq(; s) ds 



< 



RCn „ 

< -E 

~ N 



Uo 



(h Ntr (.,s), \d N)r (.,s)\ 2 )ds 



< CN^E 



[ \\dN,r(; 
JO 



s)\\ 2 ds 
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Using Doob's inequality in the form (34) we get for all T £ [0, T*]: 



E 



E 



R 



sup |M^, r (t)| 
7^te[0,T] 

' R 

^p _\Mft r (t)\ 



ds 



r=l 



*e[o,r] 



r=l 



sup ||rfiv,T-(-, 
te[o,T] 



(35) 



With (9) the analogous estimate of Mjy yields the same result as for M^ a except for the constant 
C. The correction terms M^ r and MjyV are much easier to handle and with (CI), (C2), (33) and 
(9) we get 



E 



R 

£ 



sup |M$ r (t)|+ sup |M&»I 
,te[0,T] te[0,T] 



< ctn?- 2 



(36) 



V. Discussion 

A. Connection to Astrophysics 

In the astrophysical context studied in [6, 14] the particles are dust grains in a star forming cloud 
of hydrogen gas. Depending on the grain size distribution these particles determine the opacity of 
the gas cloud and influence thermodynamic properties by emission and absorption of heat radia- 
tion. They also interact with the gas via hydrodynamic drag and influence chemical reactions via 
catalysis. A realistic expression for the collision rate of two ball shaped grains with radii l r and l q is 
a rq (x, t) ~ (l r + l q ) d ~ 1 g(\v r (x, t) — v q (x, t)\), proportional to the cross section and depending on 
the relative velocity of the two particles (cf. [14]). The material coefficient e rq i should also depend 
on the relative velocity of the two collision partners. If the latter is high, shattering collisions with 
several outgoing particles of masses mi < m r , m q are prefered. If it is low, the colliding particles 
are more likely to coagulate, so that there is one outgoing particle of mass mi = m r + m q . For a 
precise form of this function from empirical data for different grain materials we refer to [6] (table 
1) and references therein. 

Of special interest in [6, 14] is the situation after a shock with very high gas temperatures and parti- 
cle densities. In this regime shattering is the dominant process, justifying (C4), and it seems natural 
to assume that all particles are Brownian (C7). The stochasticity in the interaction coming from 
regularity conditions on W\ in (M3) also appears to be reasonable. The effective field interaction 
introduced in (M4) is a simplification we have to make in order to include the space dependence 
in a rq , which cannot be included in a model with pair interactions (see sect. V.B). The cut-off of 
the interaction rate in (M4) prevents a divergence due to high particle concentrations and seems 
natural, as for real grains the density is limited due to the positive particle volume. As the focus 
of this paper is on the interaction and not on transport terms in (MA), we left out the complicated 
hydrodynamic drag interaction between particles and gas in (M2). We note that all realistic features 
are covered by our derivation, except for the velocity dependence of e rq i, which can be included in 
a direct generalization explained in the next subsection. 

The moderate interaction scaling is technically important, as seen in the proof, and can also be 
interpreted on a physical level. In the limit N — > oo the scaling of the mean distance between par- 
ticles is given by A^v ~ N^ 1 ^ and due to (5) the interaction radius scales like ~ N~Pl d with 
(3 < 1. It vanishes more slowly with than An does and the number of interaction partners of a 
particle diverges in the limit N — > oo. In contrast to the hydrodynamic scaling with A at ~ nv this 
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introduces a self-averaging effect and the influence of the interaction partners is determined by the 
local particle density. In the astrophysical gas cloud after a shock the cross section for a collision 
interaction can be thought of being effectively enlarged by strong Brownian motion of the particles. 
Together with high particle densities this leads to a large number of interaction partners and justifies 
the moderate scaling in our microscopic model (M3), whereas the scaling in the Boltzmann-Grad 
limit [8] is not appropriate in this regime. 

B. Modifications of the microscopic model 

As explained above a modified macroscopic model (MA'), where the material coefficients e rq i 
depend on the relative velocity of the collision partners is more realistic. As the microscopic mate- 
rial coefficients e rq i G No are integer numbers they cannot depend on (x, t), because they have to 
be Lipschitz continuous for our proof (see (C4)). So we define analogous to (6) and (8) a proper set 
of m different, possible collision events 

( a *N,r q y : = ^N,rq XM(N,r,s)(k) a l Njrg (X^(s), s) dsj , i = 1, . . . ,m, 

with rates 

j N , rq (Xk(t),t):=mm{c a , 4,(Xft(t),t)((SW,,(t) * W N )(X&(t)) - 6 r , q W N @)/N)} 

and corresponding outcomes e l rql . So the process with the most probable outcome for the relative 
velocity of the collision partners can be given the highest rate, whereas the others are small. With 
this microscopic model (MF) we introduce a dependence of the collision outcome on the relative 
velocity. To obtain compatibility with the macroscopic model (MA) certain conditions on a % rq and 
e lqi nave t0 b e satisfied, and our proof of convergence applies with some minor changes. 
Another modification of the microscopic model is to include pair interactions. Here we have to 
define a process for every pair of particles with k / j: 

°iV>,(*) := Pff.rq a N,rq( S ) ds ^j for k < j 

forfc>i, 

with rates 

«aUW := ( 1 ~ S kj)XM(N,r,t){k) XM(N,q,t)U) W N {X%{t) - X 3 N {t)) a rq {t). 

Our proof of convergence can be applied, but works only for spatially independent macroscopic 
collision rates a rq (t). It assures convergence only up to a stopping time, as long as the summed 
rates J2jeM(N rt)°^N ' rqtt) 316 DOun ded uniformly in TV by some predefined constant (cf. (M4)). 
For this reason the first modification seems to be more attractive, as it covers all realistic features 
explained in section V.A. Nevertheless pair interactions are more realistic descriptions of collision 
events and coagulation could be included in such a model. 

C. Conclusion 

In this paper we specified a microscopic particle model (MI) from which we rigorously derived 
Smoluchowski's equation in the space dependent form (MA). Using the technique of the moderate 
limit developed in [1 1], we could estimate the critical fluctuation terms with martingale inequalities. 
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This scaling assumption is a good approximiation of real systems in the regime of high tempera- 
tures and particle densities, which has been of interest in the study of interstellar gas clouds after 
shocks. 

In the framework of the moderate limit, a further interesting question is the asymptotic behaviour 
of fluctuations for large system sizes and the formulation of a central limit theorem for this prob- 
lem. One can as well try to eliminate some of the technical conditions for the proof, such as the 
restrictions to the scaling parameter in (CI) or the finite number of particle sizes. There is also hope 
to proof a convergence result for a microscopic model with pair interactions without the constraints 
mentioned in V.B, by using a suitable Sobolev-norm. 

It would be certainly of most interest to derive Smoluchowski's equation in the hydrodynamic limit, 
but this task cannot be achieved with the methods used here. Nevertheless we could prove the valid- 
ity of the spatially inhomogeneous equation in a regime, which is of great interest in astrophysics, 
and where previous derivations do not apply. 
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